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Abstract 

We  examine  the  structure  and  stability  of  an  ultrathin  ceramic  film  coating  a  metal  substrate,  specifically,  an  a- 
alumina,  A1203,  film  grown  on  fee  nickel.  This  metal-ceramic  interface  may  play  a  role  in  materials  failure  of  current 
combustion  engine  thermal  barrier  coatings  (TBC’s).  We  characterize  the  (0001)  surface  of  ot-Al203  and  study  the 
effect  of  increasing  alumina  film  thickness  on  the  alumina/nickel  interface  using  periodic  slab  density  functional  theory 
within  the  generalized  gradient  approximation.  Since  Ni  forms  stable  alloys  with  Al,  it  is  not  obvious  whether  the  bonds 
between  Ni  and  alumina  will  be  Ni-Al,  Ni-O,  or  both.  Interestingly,  our  calculations  indicate  that  the  preferred 
bonding  mode  depends  on  the  thickness  of  the  alumina  film.  Namely,  for  one  monolayer  of  alumina,  the  alumina 
appears  amorphous  and  both  Ni-O  and  Ni-Al  interactions  take  place,  while  for  two  and  three  monolayers,  Ni-O 
interactions  decrease  and  Ni-Al  bonds  become  more  pronounced.  By  studying  the  effect  of  increasing  alumina  thickness 
on  the  Ni  substrate,  we  observe  a  marked  decrease  in  the  work  of  adhesion  for  thicker  alumina  coatings.  This  provides  a 
new  atomic-scale  explanation  for  the  observed  increase  in  spallation  with  increasing  thickness  of  oxide  layer  (alumina) 
that  forms  during  preparation  and  operational  cycling  of  TBC’s.  The  thickest  alumina  layers  energetically  prefer  intra¬ 
ceramic  bonding  over  A1203 /Ni  interface  formation.  Connections  to  metal  catalyst-oxide  support  interfaces  are  also 
discussed.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Interfaces  between  metals  and  ceramics  arise  in 
numerous  practical  applications,  including  elec¬ 
tronics,  heterogeneous  catalysis,  oxidized  alloys, 
electrochemical  fuel  cells,  gas  sensors,  and  thermal 
protective  coatings.  From  a  more  fundamental 
standpoint,  these  interfaces  are  interesting  simply 
because  they  are  poorly  characterized  and  their 
intrinsic  nature  at  the  atomic  level  is  poorly  un- 
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derstood.  Experimental  techniques  to  probe  such 
interfaces  have  advanced  such  that  spatially- 
resolved,  high-resolution  characterization  of  in¬ 
ternal  interfaces  is  possible  for  select  cases.  At  the 
same  time,  increased  computational  power  allows 
simulation  of  larger,  more  realistic  systems.  Ac¬ 
cordingly,  the  characterization  of  heterogeneous 
materials  interfaces  has  increasingly  become  a  to¬ 
pic  of  intense  research  efforts  over  the  past  decade. 

A  number  of  experimental  techniques  for  mea¬ 
suring  interfacial  binding  exist  in  the  literature. 
Gupta  et  al.  developed  a  measure  of  interface 
strength  whereby  a  laser  pulse  is  converted  to  a 
pressure  pulse  sent  through  the  substrate  towards 
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the  coating  interface.  The  reflected  tensile  wave 
from  the  coating  free  surface  results  in  spallation 
(deadhesion).  The  critical  stress  amplitude  that 
accomplishes  the  spallation  can  be  determined  via 
computer  simulation.  This  technique  also  has  the 
ability  to  determine  interfacial  toughness  for  cer¬ 
tain  systems  [1,2].  Similar  spallation  techniques 
have  been  employed  by  other  groups  [3,4]  and  have 
been  used  by  Gupta  et  al.  to  measure  the  tensile 
strength  of  interfaces  at  elevated  temperatures  [5]. 
Thin  film  multiple  cracking  [6]  and  tensile  stress- 
strain  measures  of  metal-ceramic  single  fibers  [7] 
can  also  be  used  to  measure  interfacial  shear 
strengths;  however,  the  stress-strain  single  fiber 
study  determined  shear  strengths  that  were  an  or¬ 
der  of  magnitude  lower  than  those  obtained  by  the 
multiple  cracking  technique  [7].  Forms  of  micro¬ 
indentation  tests  as  a  measure  of  interfacial  shear 
strength  are  useful  [8,9],  but  inhomogeneities  may 
require  numerous  measurements  of  this  type  to 
really  obtain  a  representative  sampling  of  the  in¬ 
terface.  Conversely,  direct  pull  measures  are  also 
used  to  measure  interfacial  adhesion  in  ceramic- 
ceramic  [10]  and  metal-ceramic  [11]  systems. 
Continuous  microscratch  techniques  are  able  to 
provide  a  qualitative  measure  of  adhesion  strengths 
in  metal-ceramic  and  other  bi-material  systems 
[12,13].  Likewise,  adhesion  strengths  for  metal- 
ceramic  interfaces  can  be  measured  via  particle 
profiles  for  solid-solid  interfaces  [14]  or  the  contact 
angle  that  a  drop  of  metal  makes  with  the  substrate 
[15].  Interfacial  toughness  can  also  be  investigated 
by  four  point  bending  tests  calibrated  via  finite  ele¬ 
ment  analysis  [16]  and  by  delamination  measure¬ 
ments  of  residually  stressed  coatings  deposited  over 
the  thin  film  interface  [17].  More  recently,  micro¬ 
calorimetry  measures  of  metal  adhesion  on  oxide 
thin  film  surfaces  have  been  performed,  which 
provide  estimates  of  bulk  adhesion  energetics  [18]. 
A  study  by  Ashcroft  and  Derby  employed  five 
different  measures  of  adhesion  and  found  that  no 
test  was  able  to  provide  quantitative  adhesion 
measures,  although  all  tests  provided  the  same 
adhesion  strength  ordering  of  the  samples  [19]. 
Thus,  although  experimental  techniques  are  able  to 
provide  qualitative  measures  of  adhesion  between 
metal-ceramic  interfaces,  quantitative  results  are 
generally  difficult  or  impossible  to  obtain. 


A  plethora  of  inherent  difficulties  arise  when 
attempting  to  experimentally  measure  the  intrinsic 
strengths  of  interfaces.  The  elusiveness  and  irre- 
producibility  of  quantitative  measurements  in  this 
field  are  primarily  caused  by  the  presence  of 
imperfections  (defects,  voids,  etc.)  created  upon 
sample  preparation.  It  is  virtually  impossible  to 
create  an  ideal  interface  and  to  reproducibly  mea¬ 
sure  the  work  of  separation  (adhesion  strength) 
of  such  an  interface.  As  a  result,  we  have  begun 
a  series  of  first-principles  theoretical  studies  de¬ 
signed  to  characterize  the  adhesion  and  bonding  at 
ideal  interfaces  between  ceramics  and  metals  [20- 
22].  Concurrent  with  these  studies,  we  also  are  in¬ 
vestigating  specific  effects  of  reduction  [23]  and 
oxidation  [24]  on  materials  properties  and  failure. 
Of  course,  such  calculations  isolate  the  atomic- 
level  interactions  and  cannot  account  for  certain 
long-ranged  relaxation  effects  and  dislocation/de- 
fect  emission  mechanisms  for  energy  dissipation 
present  in  physical  fracture  measurements.  Ac¬ 
cordingly,  the  adhesion  values  presented  here  may 
generally  predict  a  “lower  bound”  on  the  interface 
adhesion. 

1.1.  Metal-ceramic  interfaces  of  thermal  barrier 
coatings 

The  thermal  barrier  coatings  (TBC’s)  used  in 
airplane  jet  engine  turbines  provide  a  practical 
example  where  optimization  of  metal-ceramic  in¬ 
terfaces  is  crucial.  Similar  materials  technology  is 
used  in  such  diverse  applications  as  coating  the  gas 
turbines  of  stationary  power  plants  to  creating 
protective  outer  coatings  for  the  space  shuttle! 
What  all  such  applications  have  in  common  is 
the  necessity  of  a  protective  coating  designed  to 
withstand  the  harsh  temperature  cycling  and  the 
potentially  oxidative/corrosive  operating  environ¬ 
ment. 

During  operation,  jet  engine  components  are 
exposed  to  harsh  cycles  of  heating  and  cooling, 
spanning  a  temperature  range  in  excess  of  1400  K. 
TBC’s  are  designed  to  protect  the  metal  superalloy 
of  an  engine  from  these  extremely  high  tempera¬ 
tures  associated  with  combustion.  An  optimal 
TBC  must  possess  a  well-matched  structure  to  that 
of  the  substrate  and  maintain  this  by  similar 
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thermal  expansion  over  the  entire  range  of  tem¬ 
peratures  accessed  during  the  thermal  cycling.  It 
must  also  be  designed  to  have  favorable  adhesion 
to  the  substrate  and  very  low  thermal  conductivity. 
In  fact,  the  thermal  conductivity  of  zirconia,  the 
favored  ceramic  for  turbine  protective  coatings,  is 
so  low  that  even  a  coating  less  that  one  millimeter 
thick  can  reduce  the  effective  temperature  to  which 
the  underlying  metal  superalloy  is  exposed  by 
several  hundred  degrees  Celsius  [25]!  In  modern  jet 
engines,  the  operating  temperature  often  exceeds 
the  melting  points  of  the  metal  alloy  components. 
Accordingly,  it  is  essential  that  the  TBC  maintain 
its  uninterrupted  protective  barrier. 

A  schematic  structure  of  a  thermal  barrier  coat¬ 
ing  is  shown  in  Fig.  1.  The  substrate  is  a  nickel- 
based  alloy  of  macroscopic  thickness.  A  common 
metal  alloy  bond  coat  for  nickel-alloy  substrates 
consists  of  Ni  with  Cr,  Al,  and  Y  and  is  generally 
~  100-1 50  pm  thick  [26].  The  bond  coat  composi¬ 
tion  may  replace  or  decrease  the  amount  of  Ni  to 
add  Fe  and  Co,  especially  in  the  case  of  steel 
substrate  applications  [27,28].  This  bond  coat 
promotes  adhesion  and  is  designed  to  serve  as  a 
gradation  layer  for  thermal  expansion  coefficient 
mismatch  between  the  top  coat  and  substrate. 
Because  the  yttria-stabilized-zirconia  top  coat  is 
essentially  transparent  to  oxygen  ion  diffusion  at 
high  temperatures,  an  unintended  thermally  grown 
oxide  (TGO)  layer,  consisting  primarily  of  A1203, 
thickens  with  repeated  thermal  cycling  [29].  Al¬ 
though  the  TGO  is  not  intentionally  fabricated  as 
part  of  the  TBC,  such  a  layer  is  anticipated  and 


provides  some  protection  against  oxidation  of  the 
underlying  alloy,  since  oxygen  does  not  readily 
diffuse  through  alumina.  A1203  generally  forms 
even  prior  to  thermal  cycling,  due  to  the  condi¬ 
tions  present  in  the  deposition  chambers  where  the 
TBC’s  are  prepared.  With  repeated  thermal  cy¬ 
cling,  the  TBC  experiences  failure  in  the  form  of 
spallation  (chipping  off  the  substrate);  the  TGO 
may  be  partially  responsible  for  this  failure  [30,31]. 

TBC  spallation  is  a  well-known  problem  that 
has  been  studied  by  a  number  of  experimental  and 
theoretical  techniques.  Quantitative  linking  of  the¬ 
ory  to  experimental  findings  has  proven  difficult. 
Several  aspects  of  TBC  failure  studies  pose  diffi¬ 
culties  from  a  theoretical  standpoint:  the  large 
length-scales  of  the  phenomenon  (e.g.,  cracking 
can  occur  on  the  micron  scale  or  larger),  numerous 
small  percentage-component  “contaminants”  or 
“dopants”  in  pure  metals  and  ceramics,  time  scales 
associated  with  failure  evolution  that  are  much 
longer  than  characteristic  atomic-level  dynamics, 
and  the  harsh  oxidative/corrosive  environment 
and  repeated  cycling  through  a  wide  temperature 
range.  To  be  computationally  feasible,  theoretical 
studies  must  purposefully  ignore  at  least  some  of 
these  complications.  Ab  initio  calculations  are 
computationally  expensive,  which  generally  limits 
their  applicability  to  systems  of  ~  100  atoms.  This 
limits  the  chemical  and  physical  complexity  of  the 
interface  system  that  can  be  explored.  However,  a 
first-principles  approach  has  the  advantage  of 
being  parameter  free,  i.e.  the  atomic  interactions 
are  “unbiased”  -  an  especially  useful  feature  for 


Yttria-Stabilized-Zirconia  Top  Coat  (-250-500  p  m) 
Thermally  Grown  Oxide  (-AljO^)  (-10-25  p  m) 

Bond  Coat  (-100-150  p  m) 

Nickel  Alloy  (macroscopic) 


Fig.  1 .  Cross-section  of  a  TBC. 
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situations  where  fundamental  interactions  are 
poorly  understood.  Accordingly,  we  perform  first- 
principles  density  functional  theory  (DFT)  calcu¬ 
lations  to  characterize  an  ideal  interface  between 
the  metal  alloy  bond  coat  and  the  TGO  to  explore 
possible  atomic-level  causes  of  spallation.  Natu¬ 
rally,  it  is  likely  that  the  real  interfaces  present  in 
the  TBC  are  considerably  more  complicated  than 
the  model  Al203/Ni  interface  presented  here. 
Since  the  exact  structure  and  composition  at  the 
interface  is  not  well-characterized  experimentally, 
we  believe  a  study  of  the  ideal  interface  between 
stoichiometric  alumina  and  nickel  is  a  useful  first- 
order  approximation  to  the  physical  interface  and 
presents  valuable  insight  into  the  atomic-scale  in¬ 
teractions  present  at  an  interface  with  such  a 
chemical  composition. 

1.2.  Choice  of  interface  matching 

At  low  temperatures,  a  variety  of  “transitional” 
forms  of  alumina  exist.  These  transitional  forms 
cannot  be  considered  true  polymorphs  of  corun¬ 
dum,  a-alumina;  heating  transformations  between 
these  forms  are  not  reversible  [32].  Upon  heating, 
thermally-formed  oxides  of  metallic  aluminum 
transform  from  y-alumina  into  the  corundum,  a- 
alumina,  structure.  This  transformation  may  occur 
from  ~1 173  to  1373  K,  well  within  the  tempera¬ 
ture  range  inside  an  operating  jet  engine.  The  co¬ 
rundum  crystals  have  been  shown  to  grow  with  the 
basal,  (0001),  plane  oriented  parallel  to  the  sur¬ 
face  for  thin  films  [33].  Since  we  are  interested  in 
the  system  where  alumina  “grows”  on  the  bond  coat 
alloy,  we  chose  to  model  the  basal  plane  in  our 
calculations.  Previous  first-principles  DFT  [34,35] 
and  semi-empirical  studies  [36]  have  shown  the 
(0  001)  plane  to  be  one  of  the  most  energetically 
favorable  surfaces  of  alumina.  These  theoretical 
predictions  and  observation  that  the  (0  0  01)  face 
of  a-Al203  grows  in  during  high  temperature  oxi¬ 
dation  of  metallic  aluminum  motivate  our  inves¬ 
tigation  of  this  surface,  despite  the  fact  that 
achieving  ideal  cleavage  along  this  direction  may 
be  difficult  or  impossible  [37]. 

In  this  paper,  we  concentrate  on  the  interface 
formed  between  alumina  and  nickel;  specifically, 
the  a-Al203(000  1)/Ni(l  1  1)  interface.  Some  de¬ 


bate  regarding  the  optimal  planes  to  match  cubic 
metals  with  the  basal  plane  of  a-alumina  is  present 
in  the  literature.  Simple  arguments  based  on  sur¬ 
face  energy  gain  with  interface  formation  led  to 
conjecture  that  the  (110)  plane  should  be  most 
favorably  joined  with  the  (0  001)  plane  of  alumina 
[38].  However,  vapor-deposition  experiments  of  Ni 
on  a  (00  01)  plane  of  A1203  indicated  the  (111)- 
(0  001)  interface  is  more  stable  [39].  This  orien¬ 
tation,  preferred  due  to  its  low  lattice  misfit,  has 
also  been  reported  in  other  studies  [14,38,40].  Our 
own  matching  tests  confirmed  that  matching  the 
Ni(l  1  1)  plane  with  a-Al203(0  00 1)  minimized 
the  interfacial  lattice  misfit  1  to  3.3%  compared 
to  10.4%  and  6.9%  misfit,  respectively,  for  the 
Ni(100)  and  Ni(l  1  0)  planes  of  comparably  sized 
interface  unit  cells. 

1.3.  Experimental  characterization  of  the  Al203/Ni 
interface 

Most  experimental  investigations  of  the  clean 
Ni/alumina  interface  have  concentrated  on  grow¬ 
ing  thin  nickel  films  on  an  alumina  substrate.  Such 
studies  are  primarily  designed  to  characterize  the 
classic  “metal-support”  interactions  of  heteroge¬ 
neous  catalysis,  specifically  for  Ni  catalysts  on 
alumina  supports.  Here,  we  study  the  inverse  sys¬ 
tem,  i.e.  alumina  on  a  nickel  substrate,  since  the 
alumina  “film”  grows  in  on  the  bond  coat  alloy 
during  thermal  barrier  coating  oxidation. 

A  commonly  used  distinction  between  types  of 
interfacial  bonding  is  the  division  into  two  cate¬ 
gories:  physical  and  chemical  interactions.  These 
interactions  for  metal-alumina  interfaces  were  dis¬ 
cussed  by  Klomp  [41],  and  later  by  Dalgleish  et  al. 
[15],  who  reported  experimental  works  of  adhesion 
for  various  transition  metal-alumina  interfaces  in 
an  article  on  the  nature  of  metal-alumina  inter¬ 
actions.  The  energetics  of  the  physical  interactions 
are  described  by  the  work  of  adhesion,  given  by 
the  Dupre  equation: 


1  We  have  previously  defined  [20,21]  an  order  parameter  fi 
that  describes  lattice  misfit,  n  =  1  —  (2Q/(A\  +^2))  where  Q  is 
the  overlap  area  and  A\  and  A2  are  the  surface  areas  of  the 
metal  and  ceramic  cells. 
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^ad  Yc  Ym  Ycm  (1) 

with  Wa d  the  work  of  adhesion,  yc  and  ym  the 
surface  free  energies  of  the  ceramic  and  metal  re¬ 
spectively,  and  ycm  the  interface  free  energy.  Wad 
can  be  measured  experimentally  via  the  contact 
angle,  0,  that  a  liquid  drop  of  metal  makes  with 
the  ceramic  surface.  The  relationship,  known  as 
the  Young-Dupre  equation,  in  this  case  is 

Kd  =  Ym(l  +COS0).  (2) 

The  Wa d  reported  for  Ni  on  a-Al203  is  827  mJ / 
m2  in  H2  at  1273  K  using  a  modified  sessile-drop 
method  whereby  TEM  investigation  of  particle 
profiles  is  designed  to  study  solid-solid  metal- 
ceramic  interfaces  [14].  This  study  found  the  Ni  to 
be  nonwetting  on  the  alumina  substrate.  The  alu¬ 
mina  surface  in  this  study  was  primarily  oriented 
along  (0  001)  planes,  but  other  planes  were  pre¬ 
sent  in  smaller  amounts.  Although  the  metal  was 
consistent  in  orienting  a  (1  1  1)  plane  parallel  to  the 
(0  001)  alumina  film,  they  did  not  observe  pref¬ 
erential  rotational  orientations  of  the  metal  parti¬ 
cles  on  the  substrate.  Another  study,  measuring 
solid-solid  interface  energies  by  particle  profile 
investigations,  determined  a  Ni/alumina  work  of 
adhesion  of  670  mJ/m2  in  vacuum  (480  and  340 
mJ/m2  in  H2  and  He,  respectively)  at  1773  K  [42]. 
In  a  more  recent  study,  Venkataraman  et  al.  [12] 
used  continuous  scratch  tests  to  find  the  work  of 
adhesion  for  a  1.0  pm  nickel  film  on  an  alumina 
(1  120)  substrate  to  be  only  110  mJ/m2  for  a  maxi¬ 
mum  substrate  temperature  of  348  K  during  de¬ 
position  of  the  thin  film.  Clearly,  the  measured 
adhesion  strengths  are  all  quite  weak,  but  they 
show  a  wide  range  of  values,  characteristic  of  the 
irreproducibility  of  such  experiments  due  to  the  im¬ 
perfect  nature  of  the  interface  formed  and  the 
range  of  methods  used  to  assess  adhesion.  In 
the  current  study,  we  examine  ideal,  well-charac¬ 
terized  interfaces  to  predict  the  intrinsic  adhe¬ 
sion  strength. 

Local  bonding  and  chemistry  at  internal  inter¬ 
faces  of  several  metal-alumina  systems  have  been 
characterized  by  spatially  resolved  electron-energy- 
loss  spectroscopy  and  high  resolution  electron 
microscopy  [43,44].  Experimental  evidence  from 
studies  of  Ni/alumina  cermets  and  Ni  deposited 


on  clean  alumina  in  UHV  shows  the  Ni/alumina 
interface  forms  direct  Ni-Al  bonds  under  reducing 
conditions  at  high  temperatures  [43,45].  The  nickel- 
aluminum  bonds  were  expected  to  form  in  tetra¬ 
hedral  units  from  EELS  and  HREM  analysis  of 
Ni/alumina  cermets  [43].  However,  the  interfacial 
ceramic  termination  may  be  sensitive  to  the  spe¬ 
cific  metal  and  sample  preparation  technique;  for 
instance,  diffusion-bonded  interfaces  (where  the 
metal-ceramic  couple  is  formed  via  thermocom¬ 
pression  of  the  two  crystals,  as  in  the  work  of  Wan 
and  Dupeux  [40])  can  differ  from  those  produced 
by  molecular  beam  epitaxy. 

In  1990,  Zhong  and  Ohuchi  found  that  Ni  de¬ 
posited  on  a  clean  a-Al203(000  1)  substrate  in 
ultrahigh  vacuum  at  room  temperature  showed  no 
sign  of  reaction  in  XPS  core  level  and  valence  band 
spectra  [45];  similar  results  using  diffusion  bonding 
techniques  had  been  previously  reported  [46].  At 
1073  K  in  UHV,  alumina  was  partially  reduced  by 
Ni  deposition  [45].  In  the  presence  of  oxygen,  a 
spinel  phase,  NiAl204,  is  formed  at  high  tempera¬ 
ture  (1073  K).  The  effect  of  this  spinel  phase  on 
metal-ceramic  interface  adhesion  is  not  entirely 
clear.  It  may  introduce  additional  stress  to  the 
system  [45],  yet  some  experimental  evidence  indi¬ 
cates  that  the  spinel  phase  could  actually  improve 
adhesion  [47]. 

1.4.  Previous  theory  investigating  the  metal-Al203 
interface 

A  variety  of  previous  theoretical  studies  have 
been  designed  to  investigate  metal-ceramic  inter¬ 
actions.  An  exhaustive  review  of  all  such  work  is 
not  possible  here.  Some  of  these  studies  are  re¬ 
viewed  in  articles  by  Finnis  [48,49],  Ernst  [50]  and 
by  us  [20].  Previous  theory  investigating  “macro¬ 
scopic”  phenomena  of  Al203/Ni  interfaces  has 
studied  the  role  of  aperiodic  interfacial  void  for¬ 
mations  leading  to  spallation  [51]  and  modeled  the 
evolution  of  stress-strain  hysteresis  and  cyclic 
plasticity  [52].  Here,  we  will  discuss  only  a  small 
subset  of  the  theoretical  work  that  has  been  per¬ 
formed  to  investigate  the  microscopic  nature  of 
metal-ceramic  systems;  specifically  we  look  at 
several  studies  related  to  metal-alumina  interac¬ 
tions. 
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An  early  study  by  Johnson  and  Pepper  [53] 
modeled  Fe,  Ni,  Cu,  and  Ag  atoms  in  contact  with 
an  alumina  surface  via  self-consistent  field  Xa 
scattered- wave  calculations  on  clusters.  These 
calculations  showed  a  chemical  bond  to  be  esta¬ 
blished  between  the  d  electrons  of  the  metal  and 
the  “nonbonding”  2 p  electrons  of  oxygen  anions 
on  the  alumina  surface.  They  concluded  that  the 
bonding  was  primarily  covalent,  with  decreasing 
strength  across  the  series  Fe,  Ni,  Cu,  and  Ag. 

Three  years  after  the  Johnson  and  Pepper  pa¬ 
per,  several  other  reports  were  published  relating 
to  the  Ni/Al203  interface.  Andersen  et  al.  used 
extended  Hiickel  tight-binding  (atom  superposi¬ 
tion  and  electron  delocalization  molecular  orbital 
theory)  cluster  models  to  study  the  bonding  of  a- 
A1203  with  Ni,  Al  and  Y  “surfaces”  and  to  look  at 
cation  vacancy  diffusion.  The  alumina  cluster  was 
similar  to  that  used  in  the  Johnson  and  Pepper 
study,  but  it  was  found  that  a  10-atom  Ni  cluster 
was  much  more  weakly  bound  than  the  single  Ni 
atom  had  been  [54].  To  study  the  Ni/alumina  in¬ 
terface,  Zdetsis  and  Kunz  performed  unrestricted 
Hartree-Fock  calculations  on  a  cluster  embedded 
in  a  large,  self-consistently  determined  point  charge 
array.  Using  only  one  explicit  Ni  atom  in  their 
Ni/alumina  “interface”  calculations,  they  found 
an  average  interface  bond  strength  of  5.3  eV  and 
charge  transfer  of  1.5  e-  from  the  Ni  to  anti¬ 
bonding  states  of  the  alumina  surface  [55].  Raatz 
and  Salahub  performed  local  spin  density  scat¬ 
tered-wave  calculations  on  clusters  designed  to 
simulate  a  Ni  cluster  on  a  y-alumina  supporting 
medium  [56].  The  structural  model  for  the  y- 
alumina  support  consisted  of  A1508  or  A1508(H) 
clusters  with  Al-O  bond  distances  set  to  1.97  A, 
the  observed  bond  length  in  y-alumina  [57].  The 
theoretical  approach  in  this  study  was  comparable 
to  the  earlier  Johnson  and  Pepper  study  except 
that  a  more  accurate  electron  exchange-correlation 
potential  was  used.  They  found  that  the  alumina 
support  served  to  decrease  the  Ni  magnetic  mo¬ 
ments  and  reported  strong  hybridization  of  the  Ni- 
3 d  0-2 p  levels.  Most  of  their  paper  was  then  devoted 
to  exploring  the  effect  that  this  had  on  subsequent 
chemisorption  of  CO. 

A  more  recent  DFT  study  modeled  a  similar 
Nion-alumina  support  and  CO  chemisorption  sce¬ 


nario  [58].  Linear  combination  of  Gaussian-type 
orbital  density  functional  calculations  using  the 
Xa  approximation  to  the  exchange-correlation 
potential  were  employed.  The  results  of  this  study 
indicated  a  very  large,  8.1  eV,  adsorption  energy 
resulting  from  a  strong  ionic  bond  formed  between 
a  single  Ni  atom  and  the  oxygen  from  the  alumina 
cluster.  For  a  Ni  cluster  of  six  atoms,  that  ad¬ 
sorption  energy  decreased  to  3  eV  per  atom  due  to 
Ni-Ni  bonding. 

In  addition  to  Ni/Al203  studies,  several  DFT 
studies  have  been  reported  on  the  Nb/Al203- 
(0  001)  interface  [59-61].  Likewise,  density  func¬ 
tional  studies  of  A1203  with  Pt  and  Ag  adsorption 
[62]  as  well  as  thin  A1203  films  on  Al(l  11), 
Ru(00  01),  and  Mo(l  10)  substrates  with  subse¬ 
quent  Pt  (and  a  variety  of  other  metals  [63])  ad¬ 
sorption  [64-66]  have  been  investigated.  Jennison 
et  al.  found  that,  for  such  metal  substrates,  thin 
alumina  films  (2-3  layers  of  oxygen)  may  prefer  a 
structure  that  does  not  correspond  to  a  conven¬ 
tional  phase  of  alumina.  Specifically,  these  films 
have  Al  ions  occupying  distorted  tetrahedral  po¬ 
sitions  (or  a  mix  of  tetrahedral  and  octahedral 
positions  similar  to  k-A1203  as  found  with  a  larger 
periodic  unit  cell)  and  have  chemisorbed  oxygen  at 
the  interface  [64,66].  Although  they  did  not  spe¬ 
cifically  calculate  any  interfaces  with  Ni  or  nickel- 
aluminum  alloys,  they  anticipated  their  results  for 
alumina  on  aluminum  would  apply  to  NiAl-type 
substrates  since  they  expect  no  Ni  to  be  incorpo¬ 
rated  into  the  alumina  film  or  to  be  in  contact  with 
the  oxygen  [64];  and  they  found  local  density  ap¬ 
proximation  (LDA)  binding  of  a  Pd  overlayer  on 
the  alumina/ Al(l  11)  to  be  similar  to  an  experi¬ 
mental  estimate  for  such  a  film  on  NiAl(l  1  0)  [66]. 
By  contrast  to  the  suggestion  that  multiple  layers 
of  aluminum  metal  may  reside  at  the  alumina/ 
NiAl  interface,  Lozovoi  et  al.’s  DFT  simulation  of 
NiAl  oxidation  shows  an  Ni-rich  surface  at  the 
interface  with  aluminum  oxide  to  be  energetically 
the  most  favorable  interface  of  the  possibilities 
they  investigated  [67]. 

Previous  ab  initio  work  investigating  Ni/alu¬ 
mina  has  been  limited  to  small  clusters  of  alumina 
with  one  or  more  Ni  atoms  [53-56].  As  seen  in 
these  previous  works,  the  Ni-alumina  bonding 
strength  going  from  a  single  Ni  atom  on  alumina 
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to  a  Ni  cluster  on  alumina  shows  a  marked  de¬ 
crease.  It  seems  reasonable  to  expect  that  this 
trend  of  decreasing  Ni-alumina  interactions  with 
increasing  Ni  cluster  size  would  hold  true  -  and  be 
most  pronounced  -  in  the  case  of  an  “infinite” 
interfacial  plane.  Our  approach  represents  such  an 
Al203/Ni  interfacial  plane,  as  a  result  of  the  pe¬ 
riodic  boundary  conditions.  In  contrast  to  previ¬ 
ous  studies,  we  chose  to  model  a  nickel  substrate 
with  increasing  alumina  “coverage”,  designed  to 
model  the  growth  of  the  TGO  layer  on  the  bond 
coat  substrate  in  TBC’s;  to  our  knowledge,  such  a 
model  has  never  before  been  studied.  Further¬ 
more,  our  calculations  include  full  relaxation  of 
the  surfaces  and  interfaces  of  Al2O3(0  0  0  1)  that,  as 
we  show  later,  has  a  large  impact  on  surface  and 
interface  structure  and  energetics. 

2.  Computational  methods 

We  performed  first-principles  spin-polarized  DFT 
[68,69]  calculations  using  the  Vienna  Ab  Initio 
Simulation  Package  (VASP)  [70-72].  These  calcu¬ 
lations  employed  periodic  boundary  conditions, 
expanded  the  valence  electron  density  in  a  plane- 
wave  basis,  and  replaced  the  core  electrons  with 
ultrasoft  pseudopotentials  [73,74].  Nonlinear  par¬ 
tial  core  corrections  to  exchange  and  correlation 
were  included  for  Ni  and  Al.  For  purposes  of 
comparison,  we  performed  both  LDA  and  gener¬ 
alized  gradient  approximation  (GGA)  calculations 
for  bulk  and  (00  01)  surface  oe-Al203  calculations. 
The  interface  calculations  were  all  performed  using 
the  GGA  (PW91)  [75,76].  Dipole  corrections  per¬ 
pendicular  to  the  interface  were  calculated  for  all 
surface  and  interface  structures.  To  compare  dif¬ 
ferences  in  interface  adhesion  between  the  LDA 
and  GGA,  we  also  calculated  adhesion  energies 
within  the  LDA  for  our  most  energetically  favor¬ 
able  interface  structures.  The  lattice  vectors  used 
for  these  calculations  were  the  same  as  those  used 
in  all  interface  calculations,  but  the  ionic  coordi¬ 
nates  were  optimized  within  the  LDA  using  the 
lowest-energy  GGA  ionic  coordinates  as  the  initial 
configuration. 

The  pseudopotential  used  in  this  study  are  those 
contained  in  the  VASP  version  4.4  database.  These 


nonlocal  pseudopotentials  are  of  the  separable 
Kleinman-Bylander  [77]  form  generated  using  the 
RRKJ  scheme  [78].  The  local  portion  of  the 
pseudopotential  is  the  all-electron  potential  that 
has  been  unscreened  with  respect  to  the  valence 
electrons  beyond  a  given  radius  and  smoothly 
matched  at  this  radius  to  a  zeroth-order  Bessel 
function  except  in  the  case  of  Al,  where  the  local 
part  corresponds  to  the  d  pseudopotential.  The 
database  contains  similar  sets,  i.e.,  both  sets  use 
the  same  parameters  described  below,  that  employ 
either  the  LSD  A  of  Perdew-Zunger  [79]  or  the 
GGA  (PW91)  [75,76]  parameterization  of  the  ex¬ 
change-correlation  functional. 

The  employed  pseudopotentials  were  generated 
in  the  neutral  Asx3d9 ,  3s23pl,  and  2s22p4  configu¬ 
rations  for  Ni,  Al  and  O,  respectively.  The  local 
component  of  the  Ni  pseudopotential  had  a 
pseudization  radius  of  0.88  A.  The  ultrasoft  Ni 
pseudopotential  was  augmented  with  two  projec¬ 
tors  for  the  s,  p ,  and  d  channels.  The  s  part  had 
both  inner  (for  the  norm-conserving  part  of  the 
augmentation  functions)  and  outer  cutoff  radii  set 
to  1.14  A.  The  p  and  d  channels  for  Ni  had  inner 
radii  of  1.14  and  1.04  A,  respectively,  and  both  p 
and  d  had  outer  pseudization  radii  of  1 .29  A.  The 
Al  pseudopotential  had  two  projectors  for  the  s 
and  p  but  only  one  for  the  d.  The  Al  local  com¬ 
ponent  corresponded  to  the  norm-conserving  d 
pseudopotential  with  a  cutoff  radius  of  1 .40  A.  The 
ultrasoft  Al  s  and  p  channels  had  inner  cutoff  radii 
of  1.25  A  and  outer  pseudization  radii  of  1.40  A. 
The  O  pseudopotential  was  augmented  with  two  s 
and  two  p  projectors  with  the  local  part  having  a 
radius  of  0.54  A.  The  O  ultrasoft  s  channels  had 
inner  radii  of  0.66  A  and  cutoff  radii  of  0.80  A;  the 
p  channels  had  corresponding  radii  of  0.82  and 
1.01  A,  respectively.  Partial  core  corrections  were 
included  for  Ni  and  Al  using  truncation  radii  of 
0.51  A  for  Ni  and  0.89  A  for  Al. 

The  calculations  consisted  of  three  primary 
steps:  bulk  alumina  and  nickel  calculations,  alu¬ 
mina^  0  01)  and  Ni(l  1  1)  surface  calculations, 
and  alumina(0  00  1)/Ni(l  1  1)  interface  calcula¬ 
tions.  For  the  bulk  calculations,  we  tested  for 
convergence  of  the  k-point  sampling  density  and 
kinetic  energy  cutoff,  and  relaxed  the  atomic  and 
cell  coordinates.  We  initially  relaxed  the  surface 
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and  interface  structures  with  a  40%  lower  kinetic 
energy  cutoff  for  the  planewave  basis.  Once  we 
obtained  the  optimized  ionic  coordinates  with  this 
smaller  basis,  we  further  relaxed  the  ionic  coordi¬ 
nates  with  the  converged  planewave  basis.  This 
allowed  faster  convergence  of  the  structures,  and 
several  test  cases  showed  this  did  not  affect  the 
final  geometries.  For  the  surface  calculations,  we 
converged  the  energies  with  respect  to  slab  and 
vacuum  thickness  and  atomic  relaxations. 

The  alumina(000 1)  surface  presents  several 
cleavage  plane  options.  In  our  calculations,  we 
imposed  the  constraint  of  having  stoichiometric 
alumina  slabs  with  equivalent  faces.  This  allowed 
two  possibilities:  a  single  Al  layer  termination,  or 
an  O  layer  termination  with  half  of  the  O  atoms 
removed  on  either  side.  Fig.  2  shows  the  hexagonal 
A1203  cell  with  an  Al-terminated  (0  001)  cut  that 
can  result  in  surface  slabs  that  maintain  the  A1203 
stoichiometry  and  present  equivalent  alumina  sur¬ 
faces  on  each  side. 

In  contrast  to  the  Al203(0  001)  surface,  the 
Ni(l  1  1)  surface  presents  a  unique  cleavage  plane. 
Since  we  are  modeling  alumina  growth  on  a  nickel 


Fig.  2.  Hexagonal  unit  cell  of  ot-alumina  showing  preferred 
basal  plane  cut  in  between  the  two  Al  atoms.  The  black  spheres 
represent  aluminum  atoms  and  the  gray  ones  oxygen  atoms. 


surface,  we  chose  to  impose  the  equilibrium  bulk 
Ni  lattice  constants  on  our  interface  unit  cells.  In 
the  ionic  relaxation  steps,  the  coordinates  of  the 
nickel  atoms  in  the  bottom  layer  of  the  nickel  slab 
were  initially  fixed  to  bulk  values.  In  the  final  re¬ 
laxation  steps,  this  restriction  was  lifted,  although 
the  resulting  ionic  relaxations  in  this  nickel  layer 
were  very  minor.  As  mentioned  previously,  the 
interface  lattice  mismatch  between  the  Ni(l  11) 
and  A1203(000  1)  in  our  calculations  was  low,  at 
only  3%  for  a  hexagonal  cell  with  a  surface  area 
of  21.4  A2.  We  initially  calculated  two  different 
translational  orientations  of  the  alumina  surface  at 
the  interface  for  each  alumina  thickness.  The  first 
was  simply  the  interface  formed  by  matching  the 
(0  001)  alumina  plane  (in-plane  lattice  vectors 
[010]  and  [T 0 0])  with  Ni(l  1  1)  (in-plane  lattice 
vectors  [Oil]  and  [1 T 0]);  the  second,  more  ener¬ 
getically  favorable  translation  was  achieved  by 
uniformly  translating  (in  the  xy  plane)  the  alumina 
film  in  this  original  structure  —1.2  A  parallel  to  the 
interface.  This  resulted  in  structures  where  a  Ni 
atom  occupied  the  position  that  the  next  Al  layer 
would  have  occupied  in  the  bulk  alumina  struc¬ 
ture.  The  two  translations  corresponded  to  initial 
positioning  of  the  interfacial  Al  ions  above  an  hep 
hollow  site  in  one  translation  and  an  fee  hollow  in 
the  second  translation.  An  energetically  unfavor¬ 
able  starting  configuration,  where  the  interfacial 
aluminum  was  positioned  directly  over  an  inter¬ 
facial  nickel  site,  was  also  investigated  for  each  of 
the  three  thicknesses  of  alumina;  however,  these 
relaxed  to  the  same  structures  as  those  obtained 
with  the  energetically  favorable  fee  hollow  trans¬ 
lation  described  earlier.  To  further  test  whether 
our  relaxed  structures  were  trapped  in  local  min¬ 
ima,  the  more  energetically  favorable  (fee  hollow) 
translation  relaxed  structures  were  annealed,  using 
Nose’s  algorithm  for  constant  temperature  DFT- 
MD  [80,81],  for  0.4-0. 7  ps  at  1300  K  using  a  time 
step  of  40.0  fs  and  a  Nose  frequency  of  11.0  THz. 
The  annealed  structures  were  then  quenched  using 
the  conjugate-gradient  algorithm  also  employed 
in  the  relaxations  of  the  initial  structures  previ¬ 
ously  described. 

Providing  a  good  means  for  setting  the  partial 
band  occupancies  for  the  metal  is  important  for 
obtaining  accurate  descriptions  of  the  total  energy. 
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We  used  two  methods  in  our  calculations.  The  first 
of  these  is  the  tetrahedron  method  with  corrections 
introduced  by  Blochl  et  al.  [82],  The  second 
method  employed  in  our  calculations  was  first- 
order  Methfessel-Paxton  smearing  [83].  For  these 
calculations,  we  chose  a  smearing  width  of  0.1  eV, 
which  resulted  in  very  small  entropy  terms  (<0.5 
meV/atom).  Both  methods  gave  the  same  total 
energies  to  within  less  than  0.5  meV/atom  for  our 
calculations. 

For  analysis  of  local  charge  densities,  projection 
spheres  around  each  atom  must  be  defined.  We 
defined  spheres  with  radii  of  1.015,  1.345,  and  1.23 
A  for  Al,  O,  and  Ni,  respectively.  These  spheres 
have  no  impact  on  the  total  energy  or  density 
convergence  in  the  calculation  but  merely  deter¬ 
mine  the  amount  of  charge  density  enclosed  within 
the  predetermined  sphere  around  each  atom  for 
analysis  purposes.  The  Al  and  O  radii  were  chosen 
to  capture  most  of  the  cell  volume  of  bulk  alumina 
with  minimal  overlap;  these  radii  described  a 
crystal  of  ~80%  ionic  character  for  the  “bulk 
structure”  expanded  by  3%  to  fit  the  nickel  sub¬ 
strate.  Estimates  ranging  from  ~60%  to  ^90% 
ionic  character  have  been  reported  previously  [84- 
87].  The  sphere  radii  for  Al  and  O  were  able  to 
capture  most  of  the  electron  density;  however,  the 
Ni  radius  resulted  in  ~1.25  e-/atom  in  the  inter¬ 
stitial  regions.  This  value  was  maintained  to  mini¬ 
mize  sphere  overlap  in  our  analysis  of  charge 
differences  between  atoms  in  the  isolated  slabs  and 
the  interfacial  structures.  Given  the  partial  occu¬ 
pancy  of  the  Ni  4s  band,  which  is  expected  to  be 
delocalized  as  it  is  in  all  metals,  it  is  not  surprising 
to  find  a  large  amount  of  the  electron  density  in 
the  interstitial  regions. 


3.  Results  and  discussion 


bulk  cases  were  then  used  for  the  surface  and  in¬ 
terface  calculations.  For  the  Ni  bulk  primitive  cell 
calculations,  we  found  a  Monkhorst-Pack  grid  of 
lOx  lOx  10  k-points  to  be  sufficient.  As  expected, 
the  A1203  calculations  were  converged  with  a 
much  lower  k-point  sampling  of  3  x  3  x  3  k-points 
for  the  primitive  cell  (3  x  3  x  1  for  the  hexagonal 
cell).  A  kinetic  energy  cutoff  of  340  eV  resulted  in 
absolute  energy  convergence  of  at  least  0.01  eV  for 
our  calculations;  the  use  of  ultrasoft  pseudopo¬ 
tentials  permitted  converged  energy  values  with 
this  relatively  small  basis  (in  fact,  Ni  was  well 
converged  with  a  lower  cutoff  of  ^300  eV).  The 
kinetic  energy  cutoff  for  the  augmentation  charge 
grid  was  chosen  to  be  560  eV.  This  value  had  been 
previously  tested  for  adequacy  with  the  O 
pseudopotential  employed  and  was  higher  than 
that  required  by  the  Al  and  Ni  pseudopotentials 
[70-72];  our  tests  using  a  significantly  higher  cutoff 
(800  eV)  for  the  augmentation  charge  basis  showed 
that  the  lower  value  (560  eV)  was  indeed  con¬ 
verged.  The  bulk  structures  were  then  optimized, 
by  nonlinear  fit  to  an  equation  of  state,  using  the 
converged  k-point  and  energy  cutoff  values.  This 
was  accomplished  by  performing  single  point  cal¬ 
culations  where  we  uniformly  scaled  the  lattice 
vectors  of  the  primitive  cells  for  Ni  and  A1203  to 
sample  the  region  within  ~5%  of  the  equilibrium 
value.  The  Ni  volume  within  DFT-GGA  was 
within  1%  of  experiment.  An  earlier  study  by 
Kresse  et  al.  reported  a  bulk  modulus  [88],  ob¬ 
tained  using  the  Ni  ultrasoft  pseudopotential  in 
VASP,  that  was  within  5%  of  the  experimental 
value  [89].  For  our  alumina  calculations,  a  non¬ 
linear  fit  to  Murnaghan’s  equation  of  state  [90] 


E(V)  = 


(Vo/Vf* 

B'o~  1 


+  1 


+  C 


(3) 


3.1.  Bulk  Al2  O 3  and  Ni 

We  tested  for  convergence  with  respect  to  k-point 
sampling  density,  plane-wave  basis,  and  equili¬ 
brium  lattice  volume,  for  bulk  Ni  and  A1203.  The 
most  stringent  of  these  converged  values  for  en¬ 
ergy  cutoff  and  k-point  sampling  density  in  the  two 


(with  E(V)  the  total  energy,  V  the  volume,  Vo  the 
equilibrium  volume,  B0  the  bulk  modulus  with 
pressure  derivative  B'0 ,  and  C  a  constant)  showed 
that  both  LDA  and  GGA  results  had  bulk  moduli 
with  <10%  deviation  from  experiment  and  equi¬ 
librium  volumes  within  3%  of  the  experimentally 
determined  value  [91]. 
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Table  1 

Comparison  of  LDA  and  GGA  unrelaxed,  Eu nrei,  and  relaxed,  Erd,  surface  energies  (mJ/m2)  and  vertical  relaxation  (A z  in  A)  of  surface 
Al  (inward)  for  three-,  four-,  and  five-layer  alumina  slabs.  Overbinding  of  LDA  results  in  higher  surface  energies  for  those  structures. 
Convergence  of  surface  energy  and  structural  relaxation  is  fairly  good  by  the  four-layer  structure  for  both  LDA  and  GGA  calculations 


a-Al203(000 1) 

LDA 

GGA 

FUnrel 

Frel 

Az 

FUnrel 

Frel 

Az 

Three  layer 

4016 

1866 

0.81 

3509 

1401 

0.84 

Four  layer 

4014 

1961 

0.73 

3509 

1496 

0.74 

Five  layer 

4006 

1957 

0.71 

3505 

1494 

0.73 

3.2.  Al203(000 1 )  and  Ni( 1 1 1 )  surfaces 

The  x,  y,  and  z  coordinates  of  all  atoms  in  the 
surface  calculations  were  free  to  relax  within  the 
fixed  lattice  vectors  determined  from  equilibrium 
bulk  calculations  and  with  inclusion  of  at  least  10  A 
vacuum.  Such  a  model  produces  two  free  surfaces, 
hence  the  surface  energy  is  calculated  by  subtract¬ 
ing  the  surface  cell  energy  from  the  bulk  energy  (for 
the  same  number  of  atoms)  and  normalizing  by 
twice  the  area  of  the  xy  plane  of  the  surface  cell.  We 
calculated  the  relaxed  surface  energy  for  Ni(l  11) 
for  3-5  layer  slabs  and  found  the  surface  energy  of 
Ni(l  1  1)  to  be  converged  within  2%  of  the  five- 
layer  value  for  even  the  three-layer  case.  For  our 
alumina  surface  calculations,  we  investigated  both 
aluminum  and  oxygen-terminated  surfaces.  After 
some  preliminary  calculations  with  the  O-termi- 
nated  surface,  we  determined  this  surface  to  be 
much  less  energetically  favorable  (2920  mJ/m2  for 
a  five-layer  slab  compared  to  1490  mJ/m2  for 
a  five-layer  2  Al-terminated  slab  within  DFT- 
GGA(PW91)  calculations)  and  focused  on  the 
stoichiometric,  Al-terminated  cuts  for  all  subse¬ 
quent  surface  and  interface  calculations.  Our 
findings  of  the  Al-terminated  Al2O3(0  00  1)  surface 
as  the  most  stable  is  consistent  with  similar  pre¬ 
dictions  from  DFT  embedded  cluster  calcula¬ 
tions  [92].  The  alumina(000 1)  relaxed  surface 
energies  were  well  converged  for  four-  and  five- 
layer  slabs.  For  three-,  four-,  and  five-layer  alu¬ 
mina  slabs,  we  found  fairly  similar  (large)  inward 


2  By  “layer”  we  mean  an  ALO3  unit  that  consists  of  a  plane 
with  three  oxygen  atoms  and  an  aluminum  on  either  side  of  that 
plane,  i.e.  AI-O3-AI.  Six  such  layers  make  up  the  hexagonal 
unit  cell  of  0C-AI2O3. 


relaxation  of  the  outermost  Al  atoms,  which  re¬ 
sulted  in  a  very  large  energy  lowering  compared  to 
the  unrelaxed  (bulk-terminated)  surfaces.  The 
surface  energies  and  vertical  inward  relaxation  of 
the  surface  aluminum  for  Al203(0  001)  three-, 
four-,  and  five-layer  slabs  are  shown  in  Table  1.  In 
Table  2  we  compare  our  LDA  and  GGA  five-layer 
slab  surface  energy  predictions  with  several  previ¬ 
ously  predicted  Al2O3(0  0  0  1)  surface  energy  values 
for  slabs  with  atomic  coordinates  both  unrelaxed 


Table  2 


Comparison  of  predicted  energies  (mJ/m2)  for  unrelaxed  and 
relaxed  surfaces  for  the  basal  plane  of  a-Al203  (reported  en¬ 
ergies  rounded  to  three  significant  figures  for  consistency) 


Unrelaxed 

Relaxed 

Empirical  potentials 

Electrostatic  and  shell  modela 

6530 

2970 

Two  and  three-body  potentials15 

5040 

2040 

Two  body  and  shell  model0 

5950 

2030 

First  principles 

HF/four  layerd 

5320 

3920 

LD A/three  layer6 

3770 

1760 

LD A/five  layerf 

4010 

1960 

GGA/five  layer8 

3510 

1490 

Experimental: 

~925h 

a  Tasker  [95]  using  empirical  electrostatic  shell-model  po¬ 
tentials  [116,117]. 

bBlonski  and  Garofalini  [96]  using  a  two-  and  three-body 
empirical  potential. 

c  Mackrodt  et  al.  [36]  using  an  empirical  two-body  potential 
[118]  and  the  same  shell-model  employed  in  Ref.  [95]  for  elec¬ 
tron  polarization. 

d  Causa  et  al.  [94],  four-layer  slab  using  2D  periodic  Hartree- 
Fock  calculations. 

eManassidis  et  al.  [34,35]  using  LDA  for  a  three-layer  slab. 
fThis  paper:  LDA,  five-layer  slab. 
gThis  paper:  GGA,  five-layer  slab. 
h  Experimental  value  extrapolated  to  0  K. 
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and  relaxed.  Although  there  is  considerable  varia¬ 
tion  in  predicted  surface  energies,  the  importance 
of  allowing  atomic  relaxations  is  evident  in  each 
study.  A  recent  DFT  study  by  Wander  et  al.  found 
that  the  specific  choice  of  exchange-correlation 
functional  has  some  effect  on  the  predicted  surface 
structure  [93].  Since  alumina  is  a  highly  ionic  ma¬ 
terial,  the  bulk-like  surface  termination  is  un¬ 
favorable  since  it  leaves  “bare”  surface  Al~3+  ions; 
thus  a  large  inward  relaxation  of  surface  aluminum 
cations  results  in  a  more  energetically  favorable 
structure  whereby  the  oxygen  anions  can  screen 
more  effectively  the  positive  charge  of  the  alumi¬ 
num  ions.  We  found  large  inward  relaxation  of  the 
surface  aluminum  ions  in  all  cases,  as  has  been 
reported  previously  [34-36,94-96].  The  Al-O  bond 
length  from  bulk  calculations  using  GGA  was 
1.87  A;  this  decreased  to  1.67  A  for  the  surface 
Al-O  bonds  upon  relaxation.  The  decrease  in  bond 
length  was  associated  with  a  0.73  A  vertical  relax¬ 
ation  of  the  surface  aluminum  (on  each  side  of  the 
slab),  concurrent  with  horizontal  expansion  of  the 
hexagonally  packed  oxygen  atoms  in  the  outermost 
oxygen  planes.  Previous  studies  that  allowed  only 
relaxation  of  the  interlayer  spacing  could  not  in¬ 
clude  this  effect  [34,35,59].  More  recent  studies  in¬ 
vestigating  A1203(00  0  1),  using  Hartree-Fock  and 
LDA  calculations,  have  reported  geometrical  re¬ 
laxations  for  9-  and  18-layer  slabs  but  did  not 
calculate  surface  energies  [62,97].  LDA  calcula¬ 
tions  by  Verdozzi  et  al.,  for  a  surface  with  18  layers, 
showed  comparable  vertical  relaxations  to  our 
calculations  for  the  outermost  layers.  Both  calcu¬ 
lations  showed  inward  vertical  relaxations  for  the 
surface  aluminum  of  0.73  A,  an  outward  expansion 
of  the  terminating  oxygen  plane  of  0.03  A,  and  the 
subsequent  Al-Al  layers  showing  a  vertical  dis¬ 
tance  decreased  by  0.21  A,  compared  to  0.22  A  in 
our  five-layer  GGA  case  [62]. 

Experimental  characterization  of  the  energetics 
and  surface  structure  of  (00  0  1)  plane  of  oe-Al203 
is  somewhat  complicated.  In  fact,  the  surface 
structure  of  a-Al203(00  0  1)  is  sufficiently  compli¬ 
cated  and  conditions  dependent  to  warrant  its  in¬ 
clusion  in  a  1997  article  by  Chame  et  al.  entitled 
“Three  mysteries  of  surface  science”  [98].  Obtain¬ 
ing  the  ideal  cleavage  plane  (the  Al-terminated 
surface  studied  theoretically)  by  fracture  studies 


using  the  double-cantilever-cleavage  technique  is 
difficult  or  impossible  [37].  However,  as  mentioned 
earlier,  the  (0  001)  surface  can  grow  in  during 
A1203  formation  via  aluminum  oxidation  [33]  and 
calculations  have  shown  the  ideal  single-Al  termi¬ 
nated  surface  to  be  one  of  the  lowest  energy  sur¬ 
faces  of  a-Al203.  Below  ~1220  K,  it  is  difficult  to 
obtain  reliable  surface  energies  of  clean  A1203  due 
to  hydroxylation.  An  expression  for  the  experi¬ 
mental  surface  energy  of  the  alumina  (00  01) 
surface  is  892(ergs/cm2)  -  0.12  T(°C)  obtained  by 
measures  at  1363°C  and  4660°C  and  assuming  a 
constant  coefficient  in  the  temperature  range  [99]. 
As  we  mentioned  in  a  previous  article  [23],  the 
reliability  of  such  an  estimate  over  the  temperature 
range  is  not  possible  to  ascertain  with  only  two 
data  points,  and  its  validity  in  extrapolation  to 
much  lower  temperatures  might  also  be  suspect. 
Nevertheless,  assuming  that  a  constant  coefficient 
is  a  reasonable  approximation,  the  experimental 
estimate  of  925  mJ/m2  for  the  (00  01)  plane  of 
oc-A1203  at  0  K  is  well  below  any  of  the  theoreti¬ 
cally  predicted  values  shown  in  Table  2,  although 
our  DFT-GGA  predicted  surface  energy  is  closest 
to  this  value. 

As  expected,  the  nickel  surfaces  showed  very 
little  atomic  relaxation.  We  obtained  a  GGA  sur¬ 
face  energy  of  1950  mJ/m2  for  the  five-layer 
Ni(l  1  1)  surface,  compared  with  2011  mJ/m2  for 
previous  full  charge  density  linear  muffin-tin  or¬ 
bital  calculations  [100]  (an  all-electron  density 
functional-based  GGA  calculation)  and  experi¬ 
mentally  determined  values  of  2380  [101]  and  2450 
mJ/m2  [102].  Our  calculations  showing  that 
Ni(l  11)  has  a  higher  surface  energy  than  A1203- 
(0  001)  are  consistent  with  the  experimentally 
observed  nonwetting  of  Ni  on  alumina  [42],  as¬ 
suming  no  significant  interfacial  chemical  reaction 
takes  place  (which  seems  reasonable,  since  no  re¬ 
action  was  observed,  based  on  both  core  and 
valence  band  photoemission  spectra  during  Ni 
deposition  at  room  temperature)  [45].  The  rela¬ 
tively  high  surface  energy  of  Ni(l  11)  means  that 
Ni  will  prefer  to  aggregate  to  minimize  the  amount 
of  surface  exposed,  whereas  the  somewhat  lower 
(based  on  our  calculations)  surface  energy  of  alu¬ 
mina’s  basal  plane  suggests  that  A1203  will  favor 
allowing  exposed  surfaces. 
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3. 3.  The  Al203(0001) /Ni(l  1 1 )  interface 

From  our  previous  explanation  of  the  non¬ 
wetting  of  Ni  on  alumina,  it  seems  reasonable  to 
expect  that  coating  would  occur  in  the  reverse 
system,  alumina  on  Ni,  since  Al203’s  surface  en¬ 
ergy  is  less  than  that  of  Ni.  We  do  find  that  A1203 
can  coat  Ni;  nevertheless,  we  observe  weak  adhe¬ 
sion  for  the  alumina/Ni  interface.  Before  em¬ 
barking  on  the  interface  studies,  a  rigorous  check 
of  the  adequacy  of  our  planewave  basis  was  made. 
We  performed  one  calculation  of  a  Ni(l  11)  sur¬ 
face  with  one  layer  of  Al203(00  01)  at  a  energy 
cutoff  of  600  eV  for  the  planewave  basis  and  800 
eV  for  the  augmentation  charge.  We  found  the 
total  energy  of  this  system  was  converged  within 
0.7  x  10-4  eV  to  that  of  our  lower  energy  cutoff, 
where  we  were  converging  within  a  tolerance  of 
0.1  x  10-4  eV,  which  assured  us  of  the  adequacy  of 
our  basis.  Uniformly  expanding  the  bulk  alumina 
crystal  to  match  the  Ni(l  11)  lattice  vectors  results 
in  a  energy  increase  of  0.16  eV  per  bulk  A1203  unit, 
or  0.03  eV/atom.  Accordingly,  all  interface  adhe¬ 
sion  energies  are  calculated  by  subtracting  the 
energy  of  the  infinitely  separated  nickel  slab  and 
the  expanded  alumina  slab  from  the  total  energy  of 
the  interface  cell  calculation.  In  other  words,  the 
adhesion  energy  is  defined  as  the  energy  lowering 
achieved  by  interface  formation  compared  to  the 
energies  of  the  two  isolated  slabs. 

We  calculated  structures  corresponding  to  two 
different  translations  of  the  oe-alumina  layer  on  the 
Ni  substrate.  We  found  that  the  interface  struc¬ 
tures  relaxed  to  distinct  structural  local  minima. 
The  translation  corresponding  to  initial  aluminum 
positioning  over  fee  hollow  sites  of  the  Ni(l  1  1) 
surface  resulted  in  the  energetically  more  favorable 
relaxed  structures.  The  initially  hep  hollow  site 
positioning  produced  somewhat  higher  energy 
(local  minima)  structures  for  the  three  cases.  The 
calculated  works  of  adhesion  for  the  translations 
differed  by  10%  for  one  layer  of  alumina,  18%  for 
two  alumina  layers,  and  29%  for  three  alumina 
layers.  This  increasing  discrepancy  with  increasing 
number  of  layers  seemed  related  to  the  increasing 
“stability”  of  the  a-alumina  structure,  i.e.,  the 
three-layer  alumina  geometry  was  less  able  to  re¬ 
arrange  for  optimal  substrate  bonding  due  to  more 


favorable  intra-ceramic  bonding.  Our  most  ener¬ 
getically  stable  interfaces  corresponded  to  trans¬ 
lation  of  the  alumina  such  that  the  starting 
configuration  had  a  Ni  atom  occupying  the  space 
where  the  next  Al  layer  would  have  been  in  bulk 
alumina,  with  the  surface  aluminum  positioned 
over  an  fee  hollow  site  as  described  earlier.  From 
this  initial  point,  we  found  that  the  work  of  ad¬ 
hesion  first  increased  with  layer  thickness,  then 
decreased:  618  mJ/m2  for  one  layer  of  A1203  on  Ni, 
943  mJ/m2  for  two  layers  of  A1203  on  Ni,  followed 
by  a  decrease  to  456  mJ/m2  for  the  three-layer  case. 

Calculations  using  the  LDA  to  DFT  resulted  in 
similar  trends  in  adhesion  energies.  The  adhesion 
energies  calculated  within  the  LDA  were  larger 
than  those  obtained  from  the  GGA  calculations, 
which  is  to  be  expected  since  LDA  is  notorious  for 
overbinding.  Within  the  LDA,  using  the  GGA 
lattice  vectors  but  LDA-relaxed  ionic  coordinates, 
the  monolayer  of  A1203  on  Ni  had  an  adhesion 
energy  of  1255  mJ/m2,  two  layers  of  A1203  on  Ni 
resulted  in  a  1506  mJ/m2  adhesion  energy,  and  the 
adhesion  energy  of  three  layers  of  A1203  on  Ni 
decreased  to  858  mJ/m2.  To  test  that  the  slight 
expansion  (<2%)  of  the  GGA  Ni  lattice  vectors 
over  the  ideal  LDA  ones  was  not  a  significant 
factor  in  this  increased  adhesion,  calculations 
using  the  ideal  LDA  lattice  vectors  were  also  per¬ 
formed.  These  yielded  similar  (within  ~10%) 
adhesion  energies  to  those  obtained  with  the  GGA 
lattice  vectors.  Although  it  is  expected  that  the 
LDA  results  overestimate  the  adhesion  energies,  it 
is  possible  that  the  GGA  overcorrects  for  the  LDA 
overbinding  and  provides  a  low  estimate.  The 
GGA  weakening  of  adhesion  relative  to  LDA  for 
metal  adsorption  on  alumina  was  reported  by 
Bogicevic  and  Jennison  in  their  study  of  various 
metal  adsorptions  on  ultrathin  alumina  films  [63]. 

Cross-sectional  views  of  the  relaxed  structures 
of  one,  two,  and  three  layers  of  alumina  on  Ni  can 
be  seen  in  Figs.  3  (the  most  energetically  favorable) 
and  4  (a  higher  energy,  local  minimum  structure). 
Top  views  of  the  interfacial  A1203  layer  on  the 
interfacial  Ni  layer  for  both  translations  of  the 
three  different  alumina  film  thicknesses  are  shown 
in  Figs.  5  and  6.  Annealing  and  quenching  the  more 
energetically  favorable  structures  did  not  result 
in  enhanced  works  of  adhesion.  The  one-layer 
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Fig.  3.  One,  two,  and  three  layers  of  A1203(000  1)  on  Ni(l  1  1)  -  most  energetically  favorable  translation  at  the  interface.  The  black 
spheres  represent  Al  atoms  and  the  dark  gray  ones  are  O  atoms.  The  lower  three  planes  of  light  gray  spheres  are  Ni. 


Fig.  4.  One,  two,  and  three  layers  of  A1203(0  00  1)  on  Ni(l  1  1)  -  less  energetically  favorable  interfacial  translation  (local  minimum). 
The  black  spheres  represent  Al  atoms  and  the  dark  gray  ones  are  O  atoms.  The  lower  three  planes  of  light  gray  spheres  are  Ni. 


annealed  structure  was  energetically  identical  to 
the  initial  relaxed  structure,  and  the  two-  and 
three-layer  annealed  structures  actually  resulted  in 
structures  analogous  to  the  relaxed  structures  from 
the  (less  stable)  hep  hollow  site  translation.  Like¬ 
wise,  as  previously  mentioned,  the  energetically 
unfavorable  starting  configuration  with  the  inter¬ 
facial  aluminum  ions  positioned  directly  on  top  of 
interfacial  nickel  ions  relaxed  to  the  same  final 
structures  as  those  obtained  for  the  fee  hollow  site 
translation. 

A  single  layer  of  alumina  transforms  to  an 
amorphous  structure  that  breaks  the  planarity  of 


the  oxygen  layer  upon  relaxation;  the  two  and 
three  layer  films  of  alumina  basically  maintain  the 
oxygen  planes  as  in  bulk  alumina.  The  nearest- 
neighbor  distances  across  the  interface  for  Ni-O 
and  Ni3-Al  are  shown  in  Table  3.  The  bond 
lengths  shown  in  this  table  correspond  to  the 
nearest-neighbor  distances  between  unique  atomic 
sites  shown  in  Figs.  5  and  6.  For  both  translations, 
the  two-  and  three-layer  relaxed  structures  had  Al 
atoms  closest  to  the  nickel  surface,  with  a  planar 
hexagonal  array  of  O  atoms  above  the  interfacial 
Al  atoms.  The  single-layer  structures  break  the 
planarity  of  the  oxygen  layer,  resulting  in  nearest 
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Fig.  5.  Top  view  of  the  A1203  layer  at  the  interface  with  Ni  for  one,  two,  and  three  layers  of  A1203(00  0  1)  on  Ni(l  1  1)  -  most  en¬ 
ergetically  favorable  translation.  The  black  spheres  represent  Al  atoms,  the  red  ones  are  O  atoms,  and  the  purple  ones  Ni  atoms.  The 
yellow  “X”  ’s  denote  representative  atoms  nearest  to  the  Ni. 


AI2  O3  AI2  O3  AI2  O3 

monolayer  bilayer  trilayer 


Fig.  6.  Top  view  of  the  A1203  layer  at  the  interface  with  Ni  for  one,  two,  and  three  layers  of  A1203(0  00  1)  on  Ni(l  1  1)  -  less  ener¬ 
getically  favorable  translation.  The  black  spheres  represent  Al  atoms,  the  red  ones  are  O  atoms,  and  the  purple  ones  Ni  atoms.  The 
yellow  “X”  ’s  denote  representative  atoms  nearest  to  the  Ni. 


approach  of  O  atoms  rather  than  Al,  with  the  Al 
atoms  somewhat  farther  from  the  interface  in  the 
“middle”  of  the  alumina  film;  O  atoms  are  also 
on  the  surface  exposed  to  vacuum,  i.e.,  farthest 
from  the  nickel  surface.  As  can  be  seen  from  this 
table,  the  higher  energy  (local  minimum)  two-  and 
three-layer  A1203  structures  resulted  in  Al-Ni  bond 


distances  that  roughly  correspond  to  Ni3-Al  bond¬ 
ing,  i.e.  the  relaxed  Al  at  the  interface  occupied  a 
position  with  nearly  equal  distances  from  the  three 
Ni  atoms  at  the  interface  (hep  hollow  site),  but  did 
not  have  as  close  an  approach  for  Ni-O.  Experi¬ 
mental  binding  energies  and  relative  composition 
ratios  suggest  Ni3Al  intermetallic  interactions  at 
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Table  3 


Nearest-neighbor  distances  across  the  alumina/Ni  interface  for  two  translations  of  a-Al2O3(0  0  0  1)  on  Ni(l  1  1).  See  Figs.  3  and  4  to  see 
the  symmetry-broken  monolayer  film  structure.  Al]  is  the  Al  that  is  nearest  to  the  nickel  surface,  Al2  is  the  Al  that  is  farther  from  the 
nickel  surface  in  the  A1203  monolayer 


Ni-O  (A) 

O, 

o2 

o3 

Al  to  3  nearest  Ni  (A) 

Lowest  energy  films 

One  layer 

2.16 

2.17 

3.39 

Al]  2.68,  2.73,  2.85 

Al2  2.69,  2.75,  2.85 

Two  layer 

2.34 

2.46 

2.58 

2.25,  2.91,  3.02 

Three  layer 

2.37 

2.50 

2.59 

2.32,  2.98,  3.15 

Higher  energy  films 

One  layer 

2.05 

2.97 

3.43 

Al]  2.67,  2.88,  3.13 

Al2  2.91,  2.95,  3.29 

Two  layer 

2.77 

2.80 

2.88 

2.56,  2.58,  2.62 

Three  layer 

2.70 

2.71 

2.76 

2.63,  2.65,  2.70 

the  Ni/Al203  interface  [45].  The  binding  energy 
and  XPS  line  shapes  also  suggest  perhaps  NiAl 
interfacial  interactions,  but  not  Ni2Al3  or  NiAl3 
[45,103].  The  lower  energy  minima  for  the  two- 
and  three-layer  structures  resulted  in  a  pronounced 
Ni-Al  rather  than  Ni3-Al  arrangement,  and 
achieved  much  shorter  Ni-O  distances.  For  com¬ 
parison,  the  experimentally  determined  nearest- 
neighbor  distances  in  the  bulk  cubic  crystals  of 
NiO  (NaCl  structure),  NiAl  (B2  structure),  and 
Ni3Al  (Ll2  structure)  are  2.09  A  [104],  2.50  A,  and 
2.52  A,  respectively  [105,106].  Thus,  we  see  that 
none  of  the  Ni-O  bonds  at  the  interface  are  as 
short  as  those  in  bulk  NiO,  whereas  some  of  the 
Ni-Al  bonds  are  even  shorter  than  those  in  bulk 
NiAl  and  Ni3Al.  To  the  extent  that  bond  length 
correlates  with  bond  strength  and  given  the  weak 
adhesion,  it  appears  that  the  long  Ni-O  interac¬ 
tions  must  play  a  significant  role  in  the  trends  in 
adhesion.  However,  from  our  projected  atomic 
charge  density  analysis  presented  later,  it  is  clear 
that  conventional  ionic  Ni-O  bonds  are  not 
formed  in  these  structures;  rather,  there  is  some 
energetic  preference  for  “close  contact”  at  the  in¬ 
terface  despite  the  lack  of  ionic  bond  formation 
across  the  interface.  The  preference  for  ionic 
bonding  is  effectively  isolated  to  interactions  be¬ 
tween  Al  and  O. 

Fig.  7  displays  the  density  of  states  (DOS)  for 
the  separated  and  combined  interface  systems.  The 
absolute  energies  have  been  aligned  according  to 


the  Fermi  level  (set  to  0.0)  for  the  Ni  and  interface 
DOS.  Since  the  isolated  alumina  is  insulating,  we 
aligned  the  top  of  its  O  2s  band  to  the  O  2s  band  of 
the  interface  DOS  for  purposes  of  comparison. 
The  DOS  of  the  interfacial  system  appears  to  be 
very  nearly  the  sum  of  the  two  isolated  systems. 
Analysis  of  charge  densities  projected  onto  atoms 
shows  limited  charge  differences  of  the  interfacial 
atoms  from  those  of  the  isolated  systems.  Table  4 
shows  the  difference  between  the  isolated  and 
combined  systems’  projected  local  charge  densities 
for  the  A1203  interfacial  layer  and  the  Ni  atoms  at 
the  interface.  We  did  not  observe  charge  donation, 
as  reported  in  previous  cluster  studies  [55],  from 
the  Ni  to  our  stoichiometric  alumina  slab.  The 
overall  increase  in  charge  for  interfacial  atoms  was 
generally  due  to  increased  localization  of  the 
charge  density  in  the  interfacial  region  compared 
to  the  more  diffuse  free-surface  density.  The  lack 
of  new  DOS  features  and  charge  transfer  suggests, 
as  reported  from  previous  experimental  observa¬ 
tions  [45,46],  that  low-temperature  formation  of 
a-Al203/Ni  interfaces  results  in  very  limited  or 
nonexistent  interfacial  reaction. 

Figs.  8  and  9  display  valence  electron  density 
plots  for  vertical  cuts  through  one,  two  and  three 
layers  of  A1303  on  Ni.  One  of  the  more  dramatic 
effects  to  be  noted  from  these  figures  is  the  equi¬ 
valent  contour  from  the  Ni  reaching  toward  the 
interfacial  aluminum,  indicated  by  the  red  arrows, 
in  Fig.  9.  The  presence  of  Ni-Al  interactions  is 
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Ni  +  1  layer  Alumina 


Ni  +  2  layers  Alumina 


Ni  +  3  layers  Alumina 


Density  of  States 


Density  of  States 


Density  of  States 


Fig.  7.  DOS  of  isolated  slabs,  right,  and  interface  systems,  left,  for  one,  two,  and  three  layers  of  alumina  on  Ni(l  1  1). 


Fig.  8.  Valence  electron  charge  density  cut  near  interfacial  oxygen  of  one,  two,  and  three  layers  of  A1203  on  Ni.  There  are  10 
equivalent  contours  shown  (spanning  0  to  1).  The  darkest  regions  correspond  to  the  highest  charge  density  values. 


Fig.  9.  Valence  electron  charge  density  cut  near  interfacial  aluminum  of  one,  two,  and  three  layers  of  A1203  on  Ni.  There  are  10 
equivalent  contours  shown  (spanning  0  to  1).  The  darkest  regions  correspond  to  the  highest  charge  density  values.  The  red  arrows 
mark  the  density  contour  from  the  Ni  reaching  towards  the  Al  at  the  interface.  The  Al  itself  is  not  apparent  since  it  has  very  little 
associated  valence  electron  density. 


also  apparent  when  analyzing  the  projected  charge 
densities  (Table  4).  Also  noteworthy  are  the  de¬ 
creasing  equivalent  contours  enclosing  Ni-O  re¬ 
gions  with  increasing  alumina  film  thickness  in 
Fig.  8.  This  suggests  decreasing  Ni-O  interactions 
as  the  alumina  film  thickens.  The  decreasing  Ni-O 
interactions  can  also  be  seen  from  the  charge 


transfer  for  the  interfacial  atoms  in  Table  4.  In 
Table  4,  charge  differences  between  the  isolated 
slabs  and  the  interfacial  structure  are  apparent  for 
NiAl  (ascribed  to  the  case  where  a  single  Ni-Al 
bond  length  is  significantly  shorter  than  the  other 
Ni-Al  bond  lengths  shown  in  Table  3)  and  M3AI- 
like  (ascribed  to  the  situation  where  three  Ni-Al 


E.A.A.  Jarvis  et  al.  /  Surface  Science  487  (2001)  55-76 


71 


Table  4 

Charge  difference  (e“/atom)  projected  onto  the  interfacial  A1203  and  Ni.  The  difference  is  taken  between  the  projected  interfacial 
atomic  charges  and  those  of  the  isolated  systems  with  the  same  atomic  positions.  In  nearly  all  cases  shown,  this  difference  is  ^  0.  The 
atomic  labels  All,  Al2,  etc.  are  the  same  as  those  shown  in  Table  3.  The  Ni  labels  refer  to  the  unique  interfacial  nickel  positions,  with 
the  Nil  referring  to  the  nickel  with  smallest  nearest-neighbor  aluminum  distance.  Ni2  and  Ni3  are  the  interfacial  nickel  positions  next 
nearest  to  the  aluminum,  i.e.,  making  up  the  ~Ni3Al.  Ni4  is  the  fourth  unique  Ni  position  in  the  interfacial  periodic  unit  and  is  farthest 
from  the  interfacial  Al,  i.e.,  does  not  participate  in  Ni-Al  bonding 


Atom 

Lowest  energy  structure 

Higher  energy  structure 

1  layer 

2  layer 

3  layer 

1  layer 

2  layer 

3  layer 

Af 

0.01 

0.06 

0.04 

0.01 

0.04 

0.02 

ai2 

0.01 

- 

- 

0.00 

- 

- 

Ch 

0.07 

0.05 

0.03 

0.10 

0.01 

0.00 

o2 

0.07 

0.03 

0.02 

0.01 

0.00 

0.00 

o3 

0.03 

0.04 

0.02 

0.00 

0.00 

0.00 

Ni! 

0.09 

0.15 

0.13 

0.07 

0.05 

0.04 

Ni2 

0.05 

0.00 

0.00 

0.04 

0.05 

0.04 

Ni3 

0.04 

0.00 

0.00 

0.04 

0.04 

0.03 

Ni4 

0.00 

0.00 

0.00 

0.00 

-0.01 

0.00 

bond  lengths  shown  in  Table  3  are  roughly 
equivalent  to  each  other)  interactions;  the  Al-Ni 
nearest  neighbors  in  the  lowest-energy  structures 
each  show  relative  electron  density  increase  in  the 
interfacial  environment  relative  to  the  isolated 
surface  electron  density.  Similarly,  in  the  higher- 
energy  structure,  there  is  a  charge  increase  for  the 
interfacial  Ni3Al  atoms  relative  to  these  atoms’ 
isolated  surface  charges.  Although  the  interfacial 
O  atoms  do  show  some  charge  increase,  their 
nearest-neighbor  Ni  atoms  show  no  corresponding 
increase  or  decrease  for  the  lowest  energy  two-  and 
three-layer  structures.  As  can  be  seen  from  this 
table,  the  difference  in  atomic  charges  between  iso¬ 
lated  surface  and  interfacial  environments  was  very 
minor  in  all  cases.  Since  there  is  negligible  Ni-O 
charge  transfer,  it  seems  plausible  that  the  de¬ 
creasing  Ni-O  charge  density  contours  in  Fig.  8  is 
related  to  the  increasing  stability  of  the  alumina 
structure  for  increasing  alumina  film  thickness. 

Although  these  density  plots  are  somewhat  in¬ 
formative,  perhaps  a  more  useful  tool  for  charac¬ 
terization  of  bonding  is  the  electron  localization 
function  (ELF)  [107,108].  This  function  is  helpful 
in  the  classification  of  chemical  bonding  charac¬ 
teristics  by  analysis  related  to  the  Pauli  exclusion 
principle.  Mathematically,  the  ELF  takes  the  form 


ELF  (r)  = 


1 

1  +  [£>(r)/A,(r)]2 


(4) 


where  D( r)  is  the  Pauli  excess  kinetic  energy  den¬ 
sity  3  and  Dh(r)  is  the  kinetic  energy  of  the  ho¬ 
mogeneous  electron  gas  for  a  density  equal  to  the 
local  density  [109].  The  ELF  value  is  0.5  for  a 
homogeneous  electron  gas  and  can  increase  to 
~1.0  for  electrons  paired  in  a  covalent  bond  or  a 
highly  localized,  unpaired  electron  in  a  dangling 
bond  [110].  ELF  values  below  0.5  are  less  readily 
understandable.  These  can  correspond  to  regions 
of  very  low  electron  density.  Figs.  10  and  11  show 
cross-sectional  ELF  cuts  of  the  Al203/Ni  system 
for  one,  two,  and  three  layers,  for  cutting  near 
interfacial  O  atoms  and  Al  atoms  respectively.  4 
Fig.  10  shows  very  little  localization  between  the 
interfacial  O  ions  and  the  Ni  surface,  especially  for 
the  two-  and  three-layer  cases.  In  sharp  contrast 


3  D(r)  =  T  +  rcorr  —  rBose>  where  T  is  the  kinetic  energy  of 
the  noninteracting  Kohn-Sham  system  (—2 the 
kinetic  energy  “correlation”  correction  component  is  Tcorr  = 
(l/2)(h2/2 m)V2p  to  ensure  T  +  rcorr  is  positive  definite,  and 
Those  is  the  ideal  Bose  gas  kinetic  energy  ((l/4)(7j2/2m)  x 
(|Vp|2/p)). 

4  The  low  ELF  values  within  the  nickel  slab  arise  from  the 
delocalized,  metallic  electron  density  of  the  Ni.  The  ultrasoft 
pseudopotential  augmentation  density,  which  primarily  affects 
the  localized  region  around  Ni  nuclei  as  can  be  seen  from  the 
density  plots,  is  not  included,  since  this  density  is  not  contained 
in  the  Kohn-Sham  orbitals.  These  plots  display  a  qualitative 
picture  of  bonding  between  Ni-Al  and  show  strong  localization 
of  valence  electron  density  on  oxygen  anions  as  expected. 
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Fig.  10.  ELF  (of  the  ultrasoft  pseudocharge  density)  cut  near  interfacial  oxygen  of  one,  two,  and  three  layers  of  A1203  on  Ni.  There  are 
eight  equivalent  contours  to  span  a  range  from  0.0  to  our  highest  value  of  ~0.85  for  these  ELF  cuts.  The  darkest  regions  are  the  region 
of  highest  ELF  value. 
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Fig.  1 1 .  ELF  (of  the  ultrasoft  pseudocharge  density)  cut  near  interfacial  aluminum  of  one,  two,  and  three  layers  of  A1203  on  Ni.  There 
are  eight  equivalent  contours  to  span  a  range  from  0.0  to  our  highest  value  of  ~0.85  for  these  ELF  cuts.  The  darkest  regions  are  the 
regions  of  highest  ELF  value. 


to  Fig.  10,  Fig.  11  shows  a  higher  degree  of  lo¬ 
calization  between  the  interfacial  Al  and  Ni, 
most  notably  for  the  two-layer  case.  It  seems  likely 
this  electron  localization  between  the  Ni-Al  is 
largely  responsible  for  the  high  Wa d  value  for  the 
two-layer  case  compared  to  the  one  and  three  layer 
IFad5S. 

3.4.  Summary  and  implications  for  thermal  barrier 
coatings 

We  have  investigated  the  (00  01)  surface  of  a- 
alumina  and  its  interface  with  Ni(l  11).  Consis¬ 
tent  with  the  experimental  observations  of  Zhong 
and  Ohuchi  [45],  we  find  limited  interaction  be¬ 
tween  the  surfaces  at  low  temperature  in  vacuum. 
Charge  transfer  across  the  interface  is  negligible 
in  our  calculations,  a  markedly  different  predic¬ 
tion  from  those  of  previous  Ni/alumina  cluster 


calculations  [55,58].  As  mentioned  earlier,  previ¬ 
ous  extended  Hiickel  tight-binding  cluster  simu¬ 
lations  showed  a  marked  decrease  in  adhesion  of 
Ni  clusters,  compared  to  that  of  Ni  atom,  on  an 
alumina  substrate  [54].  This  is  consistent  with  the 
results  of  recent  studies  that  have  investigated 
increasing  metal  coverage  from  1/3  ML  to  1  ML 
of  non-Ni  metals  on  alumina  where  they  find 
ionic  bonds  dominate  at  low  metal  coverage  but 
that  the  attraction  is  a  result  of  polarization  at 
high  metal  coverage  [63].  We  find  unexpected  and 
unconventional  changes  in  the  interfacial  struc¬ 
tures  and  interactions  going  from  one  to  three 
layers  of  alumina  film  thickness  on  Ni.  We  find 
that  the  Ni-O  interactions  are  more  pronounced 
in  the  monolayer  structure,  Wad  is  strongest  for 
the  two -layer  case  due  to  localized  Ni-Al  inter¬ 
actions,  while  the  three-layer  structure  exhibits 
much  weaker  Al203-Ni  interactions  at  the  inter¬ 
face,  instead  preferring  intra-ceramic  interactions. 
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These  subtle  and  complex  collective  effects  would 
likely  prove  difficult  or  impossible  to  model  with 
simple  analytic  potentials.  By  looking  at  the  re¬ 
lated  system  of  increasing  alumina  thickness  on  a 
nickel  substrate  and  finding  a  dramatic  decrease 
in  adhesion  for  the  thickest  alumina  coatings 
studied,  our  results  complement  the  findings  of 
the  Ni  cluster  study  and  provide  further  insight 
into  the  nature  of  the  interfacial  interactions  for 
AI2O3/  Ni.  The  weak  adhesion  for  this  interface  is 
supportive  of  clean  Al203/Ni  displaying  weak 
metal-catalyst  support  interactions.  It  is  likely 
that  any  strong  interactions  that  exist  between  the 
support  and  active  phase  are  primarily  due  to 
intermediate  oxide  and  aluminate  formation,  as 
indirectly  suggested  by  de  Korte  et  al.  in  their 
XRD  and  reduction  behavior  measurements  of 
coprecipitated  catalysts  [111]. 

As  mentioned  earlier,  the  actual  interface  be¬ 
tween  the  bond  coat  alloy  and  alumina  in  a  TBC 
system  might  be  highly  complex.  We  believe  the 
interface  studied  here  to  be  a  reasonable  ap¬ 
proximation  to  the  interface  resulting  from  TBC 
bond  coat  oxidation  since  excess  oxygen  is 
available  to  react  with  any  aluminum  that  mi¬ 
grates  to  the  ceramic-metal  interface,  and  the 
percentage  composition  of  aluminum  in  the  bond 
coat  is  small.  Although  some  progress  has  been 
made  using  DFT  calculations  to  predict  the  oxi¬ 
dation  of  NiAl  [67]  and  a  nonstoichiometric  de¬ 
scription  of  the  alumina  might  prove  desirable 
[64],  an  exact  description  of  interfaces  in  complex 
systems  is  still  far  from  achievable  with  current 
computational  size  limitations.  Nevertheless,  we 
believe  this  in-depth  investigation  of  an  ideal 
Ni(l  11)  interface  with  0C-AI2O3  provides  valuable 
insight  into  adhesion  trends  and  bonding  in  this 
system. 

Recent  experiments  have  shown  delamination 
between  the  bond  coat  alloy  and  alumina  in  failed 
TBC’s  as  well  as  fracture  between  the  zirconia  top 
coat  and  alumina  [112-114].  Our  work  provides 
new  insight  into  the  atomic-level  cause  for  this 
bond  coat-alumina  delamination  with  increasing 
alumina  thickness.  Namely,  the  strain  from  ther¬ 
mal  cycling  and  bond  coat  oxidation  may  readily 
permit  delamination  between  the  alumina  and 
nickel  due  to  the  very  weak  interfacial  adhesion. 


Our  previous  characterization  of  the  Zr02/Al203 
interface  led  to  the  conclusion  that  this  ceramic- 
ceramic  interface  is  also  weakly  adhered  [21].  Ac¬ 
cordingly,  limiting  alumina  formation  may  be  a 
key  factor  in  extending  the  operating  lifetime  of 
the  TBC. 

Excessive  bond  coat  oxidation  leading  to  thick 
alumina  scale  formation  is  detrimental  to  TBC 
performance  lifetimes;  finding  means  to  limit  this 
process  should  help  to  solve  this  problem.  Dis¬ 
covering  chemical  means  to  protect  the  superalloy 
from  oxidation  by  forming  slow-growing  oxides 
with  favorable  adhesion  may  be  the  critical  aspect 
of  materials  technology  that  will  allow  further 
advances  in  TBC  performance.  Recent  experi¬ 
mental  work  by  Pint  et  al.  investigated  the  effect  of 
reactive  elements  such  as  Pt  on  alumina  scale  ad¬ 
hesion  [30].  They  found  that,  although  Pt  addition 
to  a  nickel  aluminide  bond  coat  did  improve  me¬ 
tal-ceramic  adhesion,  it  was  not  as  successful  as 
certain  other  elements,  i.e.,  Zr  and  Hf.  Pint  et  al. 
speculate  this  was  because  it  was  not  as  effective  in 
reducing  the  oxide  scale  growth  rate  or  changing 
the  0C-AI2O3  crystal  structure.  Nevertheless,  Pt 
addition  to  bond  coatings  in  TBC’s  has  become 
rather  popular.  Another  study  found  that  deposi¬ 
ting  a  layer  containing  Pt  between  the  zirconia  top 
coat  and  the  bond  coat  alloy  improved  oxidation 
resistance.  They  interpreted  this  improvement  as 
due  to  Pt  decomposition  of  alumina  at  high  tem¬ 
peratures  [115].  Of  course,  a  natural  approach  to 
limiting  alumina  scale  growth  would  be  to  limit  the 
Al  content,  and  hence  A1203  formation,  in  the 
bond  coat  alloy  itself.  We  are  currently  investi¬ 
gating  candidate  metals  with  which  to  replace 
some  or  all  of  the  aluminum  content  in  the  boat 
coat.  This  may  result  in  TBC’s  that  can  be  fabri¬ 
cated  with  the  same  methods  as  those  currently 
employed,  but  which  will  have  longer  performance 
lifetimes  due  to  increased  oxidation  resistance  and / 
or  improved  adhesion  of  the  zirconia  to  the  bond 
coat  alloy. 
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